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Vortex physics plays an important role in the study 
of unconventional superconductors. In this short note, 
based on Eilenberger theory we study the temperature 
(T) dependence of lower critical field H c i(T) of vortex 
states in anisotropic superconductors. 

Recent developments of experimental technique make 
us possible to observe H c i(T) exactly, and discuss it 
in the relation to the mechanism of unconventional su- 
perconductivity, such as in new iron-based superconduc- 
tors 1 ' 2 ' and URu2Si2. 3 ' In traditional Ginzburg-Landau 
(GL) theory, H c \ is given by 4 ' 
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H cl cx X(T)- 2 (\nK + c ) 



(1) 



with penetration depth A, GL parameter k, and a con- 
stant Co. In the London theory A(T) = A(T = 0) + rJA(T) 
behaves as SX(T) cx exp(— A/A^T) at low T in the s- 
wave pairing, reflecting superconducting gap A. In d- 
wave pairing with line nodes, SX(T) oc T at low T in the 
clean limit. These indicate that H c i(T) depends on the 
pairing symmetry of anisotropic superconductors. 

We note that GL theory is a phenomenological theory 
valid near the transition temperature T c . Thus it is not 
clear whether the above-discussion on H C \(T) is quan- 
titatively valid. Therefore, it is expected that H cl (T) is 
evaluated by Eilenberger theory, which is quantitatively 
reliable in vortex states even far from T c . To study con- 
tributions by the pairing symmetry, we calculate H c i(T) 
for s-wave pairing with full gap and d x 2_ y 2-wetve pair- 
ing with line nodes, as typical examples, by quantitative 
Eilenberger theory. The previous work for s-wave pair- 
ing was done in a single vortex. 5 ' Our calculation is per- 
formed in vortex lattice. We also study H C \{T) for chiral 
p±-wave pairings, to see dependences on the chirality di- 
rections, i.e., parallel or anti-parallel to applied fields. 

In this study, for simplicity, we use isotropic cylindri- 
cal Fermi surface k = fcp(cos0,sin0) and magnetic fields 
are applied to the c direction. The quasiclassical Green's 
functions g(u) n , k, r), /(w„, k, r), and p(ui n , k, r) are cal- 
culated by the Eilenberger equation 

{w n + v(V + iA)}/ = A0 9) 

{w n -v(V-iA)}/t = AV*0, (2) 

g = (1 — J/^) 1 / 2 in the vortex lattice state, with the 
selfconsistent conditions of pair potential 

A(r) =g N T £ (^(k) (/ + /+' 

0<lx> rl <tJ cu t 



(3) 



and the vector potential 



2T 

Vx(VxA) = T V (vlmg) 



(4) 



0<w„ 



in Eilenberger unit, 6,7 ' with Matsubara frequency cj n , 
(ffo^o)" 1 = InT + 2TJ2 Q<LJn ^ ut u ~\ where v - k/fc F 
is the direction of Fermi velocity vp, r is the center- 
of-mass coordinate, and (■ • • )k indicates the Fermi sur- 
face average. We use w cu t = 20£;bT c . The internal field 
B(r) = B + V x a(r) is related to the vector poten- 
tial A(r) = |B x r + a(r) in the symmetric gauge, 
where B = (0, 0, B) is a uniform flux density. The pair- 
ing function is defined as </>(k) = 1 for s-wave pairing, 
<^>(k) = v / 2cos2# for d x 2_ y 2-wave pairing. In the chi- 
ral p-wave pairing, we consider two-component order pa- 
rameter A_|_(r)0 + (k) + A_(r)0_(k) instead of A(r)0(k), 
where 0±(k) = e ±lS . 8 ' In the p + - (p_-) wave pairing, A + 
(A_) is main component with singular vortex, and A_ 
(A + ) is passive component induced around vortices. 

Our calculation is done for k = 2 and triangular vortex 
lattice. We iterate calculations of eqs. (2)- (4) under given 
B, and obtain selfconsistent vortex solutions for spatial 
structures of A(r), A(r), and quasiclassical Green's func- 
tions, as done in previous works. 7 ' 8 ' Using the solutions, 
we calculate the external magnetic field H by 



H = B 



B(r) - B) 
T ^ //l 



k B 



E 

0<u„ 



-Re 



IB 

(/tA + /A*) 9 



t Re{g - 1} 



(5) 



which is derived by Doria-Gubernatis-Rainer scal- 
ing, 7 ' 9, 10 ' and (• ■ ■ ) r indicates the spatial average. Mag- 
netic fields are in unit of Bq = 4>q/2ttRq with the flux 
quantum <pQ and Rq = hv-p /2i:k^T c . 

Figure 1(a) present magnetization curves of B as a 
function of H at some T for ri a .2_j / 2-wave pairing. There, 
H c i is defined as onset of B. In Meissner states at H < 
H c i, B = 0. In Fig. 1(b), we present H c i(T) as a function 
of T for some pairing symmetries, and we replot them as 
H cl {T)/H cl {T = 0.05T C ) in Fig. 1(c) to compare the 
T-dependence each other. 

First, we discuss the differences between the s-wave 
and the d-wave pairings. H c ± in c?-wave pairing is smaller 
than that in s-wave pairing, because the condensation 
energy of d-wave pairing is weaker due to the line node 
contributions, compared to that in the full-gap s-wave 
pairing. H c ± is related to the energy for creation of a vor- 
tex in Meissner states. 4 ' As for T-dependence, reflecting 
low energy excitations by line nodes, H c i(T) in d-wave 
pairing decreases rapidly at low T, compared with s-wave 
pairing. 

To discuss quantitative validity of the relation in eq. 
(1), in Fig. 1(c) we also present X^ ndon (T) given by 

|A0| 2 v 2 



^London X ^ E 



(U, 2 + |A^| 2 )3/2 



(6) 



in London theory, where T-dependence of A is deter- 
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Fig. 1. (Color online) (a) Magnetization curve of B as a function 
of applied field H at T/T c = 0.1, 0.3, 0.5 and 0.7 for d-wave 
pairing, (b) T-dependence of H c i(T) for s~, d-, p+- and p_- 
wave pairings, estimated by Eilenberger theory, (c) H c i(T) in 
(b) is replotted as normalized H cl (T)/H cl (T = 0.05T C ). We also 
present normalized Aj~^ ndon (T) for s- and d-wave pairings. 

mined by gap eq. (3) in uniform states. In the s-wave 
pairing, as shown in Fig. 1(c), normalized H C \(T) in 
Eilenberger theory appears smaller than A~~o ndon (T), an d 
shows decreases even at low T. This indicates that the 
vortex core energy still has T-dependence at low T, 
rather than saturation expected by X-L OU don{T). This may 
include the contribution of vortex core shrink on low- 
ering T by Kramer-Pesch effect. 11 ) On the other hand, 
in the d-wave pairing, H C \{T) in Eilenberger theory is 
higher than ALo ndon (T). Thus, T-dependence of the core 
energy is weaker than estimate by ALondon(T). This is 
an opposite effect to the s-wave pairing case, and indi- 
cates that the estimate of core creation energy is not 
simple in cZ-wave pairing because we have to consider 
both contributions inside and outside of vortex cores. 
The latter is contributions by quasiparticles extending 
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toward node-directions. 12 ) These behaviors of H C \{T) is 
also confirmed for k = 6.9. We expect that the relation 
in eq. (1) will be examined in experiments, comparing 
H cl (T) with \- 2 (T). 

Next, we study H c i(T) in chiral p-wave superconduc- 
tors. H c i(T) in p_-wave pairing is smaller than that in 
p+-wave pairing. This difference in quantitative estimate 
is consistent to previous results by phenomenological GL 
theory. 13 ' 14 ^ In chiral p- wave superconductors, opposite 
chiral component is induced around vortices of main chi- 
ral component, and core energy becomes smaller by the 
induced component. Compared with p + -wave pairing, 
the induced component is larger in p_-wave pairing, and 
the core energy is smaller, making H c \ smaller. If do- 
mains of p + -wave pairing and p_-wave pairing coexist 
at a zero-field, on increasing fields vortices penetrate at 
lower H c \ only into the p_-wave domain, where chiral- 
ity is antiparallel to the applied field. 14 -' As for the In- 
dependence, in Fig. 1(c) we see that normalized H c \ both 
for p + - and p_-wave pairings have similar T-dependence 
to that in s-wave pairing. This is reasonable, because 
p±-wave pairing with \<f>±\ = 1 has full gap, as in s-wave 
pairing. 

In summary, we quantitatively estimated different T- 
dependences of H c \ between s-wave and <i-wave pair- 
ings by Eilenberger theory. The T-dependences of H c i(T) 
show quantitative deviation from AL rad n (T). We also 
studied differences of H C \(T) between p + and p_-wave 
pairing in chiral p-wave superconductors. We expect that 
future experimental studies will confirm the relations of 
H c i(T) and the pairing symmetry in various anisotropic 
superconductors. 

We would like to thank K.M. Suzuki and K. Inoue for 
fruitful discussions, and their supports for calculations. 

1) R. Okazaki, M. Konczykowski, C. J. van der Beek, T. Kato, 
K. Hashimoto, M. Shimozawa, H. Shishido, M. Yamashita, M. 
Ishikado, H. Kito, A. Iyo, H. Eisaki, S. Shamoto, T. Shibauchi, 
and Y. Matsuda: Phys. Rev. B 79 (2009) 064520. 

2) Z. Pribulova, T. Klein, J. Kacmarcik, C. Marcenat, M. Kon- 
czykowski, S. L. Budko, M. Tillman, and P. C. Canfield: Phys. 
Rev. B 79 (2009) 020508(R). 

3) R. Okazaki, M. Shimozawa, H. Shishido, M. Konczykowski, Y. 
Haga, T. D. Matsuda, E. Yamamoto, Y. Onuki, Y. Yanase, T. 
Shibauchi, and Y. Matsuda: J. Phys. Soc. Jpn. 79 (2010) 084705. 

4) A. L. Fetter and P. C. Hohenberg: in Superconductivity, ed. R. 
D. Parks (Dekker, New York, 1969), Chap. 14. 

5) L. Kramer and W. Pesch, Solid State Commun. 12 (1973) 549. 

6) U. Klein: J. Low Temp. Phys. 69 (1987) 1. 

7) M. Ichioka and K. Machida: Phys. Rev. B 76 (2007) 064502. 

8) M. Ichioka and K. Machida: Phys. Rev. B 65 (2002) 224517. 

9) K. Watanabe, T. Kita, and M. Arai: Phys. Rev. B 71 (2005) 
144515. 

10) M.M. Doria, J.E. Gubernatis, and D. Rainer: Phys. Rev. B 41 
(1990) 6335. 

11) L. Kramer and W. Pesch: Z. Phys. 269 (1974) 59. 

12) M. Ichioka, A. Hasegawa, and K. Machida: Phys. Rev. B 59 
(1999) 184; M. Ichioka, A. Hasegawa, and K. Machida: Phys. 
Rev. B 59 (1999) 8902. 

13) R. Heeb and D.F. Agterberg: Phys. Rev. B 59 (1999) 7076. 

14) M. Ichioka, Y. Matsunaga, and K. Machida: Phys. Rev. B 71 
(2005) 172510. 



Short Note 



